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ABSTRACT 

We determine central values and radial trends in the stellar populations of the bulges of a sam- 
ple of 28 edge-on SO-Sb disk galaxies, 22 of which are boxy/peanut-shaped (and therefore 
barred). Our principal findings are the following. (1) At a given velocity dispersion, the cen- 
tral stellar populations of galaxies with boxy/peanut-shaped bulges are indistinguishable from 
those of early-type (elliptical and SO) galaxies. Either secular evolution affects stellar popu- 
lations no differently to monolithic collapse or mergers, or secular evolution is not important 
in the central regions of these galaxies, despite the fact that they are barred. (2) The radial 
metallicity gradients of boxy/peanut-shaped bulges are uncorrected with velocity dispersion 
and are, on average, shallower than those of unbarred early-type galaxies. This is qualitatively 
consistent with chemodynamical models of bar formation, in which radial inflow and outflow 
smears out pre-existing gradients. 

Key words: galaxies: abundances — galaxies: bulges — galaxies: elliptical and lenticular, cD 
— galaxies: spiral 



1 INTRODUCTION 

Among the observational clues and tests for theories of galaxy for- 
mation and evolution are the well-known empirical correlations 
between the optical absorption line strengths of early-type galax- 
ies (ellipticals and SOs, ETGs hereafter) and their dynamical or 
photometric properties (e.g. |Terlevich et al.|[l981| |Bender et al.| 
[T993] [Kuntschne r|2000| >. For increasing stellar velocity dispersion 
or mass, central and global observations of the line strengths (and 
therefore stellar populations) of ETGs show them to be older and 
more metal-rich, and more abundant in a-elements relative to Fe 
(e.g. |Thomas et a l. 2005; Kuntsc hner et al.|2 010l. Radial trends in 
stellar populations are a powerful additional discriminant. Mono- 
lithic collapse models (e.g.lEggen, Lynden-Bell, & Sandage|1962| 



Larson 1974; Carlberg 1984) predict steep negative metallicity gra- 
dients of -0.35 < A[Z/H] < -1.00 Hierarchical models predict 
that these gradients become generally shallower as they are diluted 
by mergers (e.g. |Wmte|i980"l|Di Matteo et al.|2009) . 

Observations of stellar population gradients in ETGs (e.g. 



Davies et al.||1993[ 


Mehlert et al. 2003 ; Kuntschner et al.||2006| 


Spolaor et al.|2009| 


Kuntschner et al.|2010||Roediger et al.|2011|> 



show that, on average, Alog(age/Gyr) w and A[a/Fe] w 0. How- 
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1 We follow convention by characterizing radial gradients with linear fits in 
logR, and using the notation AX = dX/dlogR, where X is [Z/H], [a/Fe] 
orlog(age/Gyr). 



ever, A[Z/H] ~ —0.2 ±0.1, and below a central velocity dispersion 
of ~ 150kms _1 (equivalent to a dynamical mass ~ 3 x 10 1() M S ), 
there is some evidence of a correlation: more massive systems 
have steeper negative metallicity gradients (e.g. Spola oret al.|2009| 
|Kuntschner et al.|2010| >. Above this characteristic mass, the corre- 
lation between gradient and velocity dispersion or mass disappears, 
but the average metallicity gradient is still negative. 

Analysis of the stellar populations of the bulges of spiral galax- 
ies is more challenging than in ETGs: bulges are fainter and em- 
bedded in disks, and they may have complicated or ongoing star 
formation histories, nebular emission, and significant dust absorp- 
tion. Nevertheless, the consensus is that the stellar populations of 
bulges in SO-Sbc disk galaxies are very similar to those in ellipti- 
cals and SOs. At a given velocity dispersion, their line strengths and 
the implied stellar populations are, on average, the same as those 
of earlier types ( |Proctor & Sansorn] |2002l|Thomas & Davies|2006{ 
Falcon-Bar roso et al.|2 006 ; MacArthur et al. 2009). In combination 
with dynamical and structural similarities, this is evidence that the 
bulges of these galaxies formed in a similar way to ETGs. How- 
ever, the good agreement between bulges and ETGs breaks down 
in the bulges of late-type spirals. For example, |Ganda et al.| ( |2007] > 
observed that the bulges of a sample of 18 Sb-Sd galaxies have 
smaller Mgb indices and larger Hj3 indices at a given stellar veloc- 
ity dispersion a . The differences may be evidence that these bulges 
formed or are currently affected by different processes to ETGs and 
the bulges of earlier-type disks, the most likely being the secular 
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and internal rearrangement of disk material (Kormen dy & Kenni-| 
|cutt|2004| >. 

The radial behaviour of the stellar populations of bulges in spi- 
ral galaxies is relatively uncertain. The major studies are Mo orthy &| 
Holtzman (2006 long-slit observations of the bulges of 38 galaxies 
ranging from SO to Sbc), |Falc6n-B arroso et al. (20061 and |Peletier| 
|et al.| (2007| S AURON integral field observations of the bulges of 
24 Sa galaxies), [Ganda et al.| < [2007| S AURON integral-field obser- 
vations of the bulges of 18 Sb-Sd galaxies), Ja blonka et al.|(2007| 
minor-axis long-slit observations of the bulges of 30 edge-on SO 
to Sc galaxies), Morelli et al. (2008 long-slit observations of the 
bulges 18 SO-Sbc galaxies) andjMacArthu r et al.| ( |2009| long-slit 
observations of the bulges of 8 Sa-Sd galaxies). Among the areas 
of agreement of these studies is the finding that most bulges have 
negative metallicity gradients. Compared to ETGs, however, there 
is rather more scatter in these gradients and more diversity in the 
structure of the radial profiles. This observation brings us on to the 
subject of bars (and, by implication, boxy/peanut-shaped bulges), 
since they may be the origin of that diversity. 

Simulations of the chemodynamical evolution of barred disk 
galaxies find that the bars drive an inflow of gas within corotation 
and outflow beyond (e.g. |Friedli et af1[l994"l |Friedli|[T998) . This 
radial transport naturally flattens pre-existing population gradients. 
Moreover, simulations by Wozniak ( 2007 ) predict local minima of 
stellar age at the ends of bars. Observational tests of these predic- 
tions are crucial, because bars are found in two-thirds of disk galax- 
ies and likely play a significant role in transforming the stellar pop- 
ulations of their hosts. Perez et al. have made the first systematic 
attempt to study radial stellar population trends in barred galaxies 
(Perez et al.|2007l [20091 [Perez & Sanchez-Blazquez|201 1| PSB11 
hereafter). They used optical spectra taken with long slits oriented 
along the bar major axes of 20 face-on and moderately-inclined 
barred SO-Sb galaxies. 

This Letter studies central values and radial trends in the stellar 
populations of bulges of barred galaxies by determining absorption 
line strengths and single stellar population (SSP) equivalent ages, 



[Z/H] and [a/Fe]. We use the Bureau & Freeman 1 1999| > sample 
of edge-on disk galaxies. The majority (22/28) of the galaxies in 
this sample host bulges that are boxy or peanut-shaped, i.e. bars 
viewed in projection (e.g. Combes et al. 1990; Kuijken & Merri- 
field|1995[|Bureau & Freeman|1999[|Chung & Bureau|2004||Kr> 



rmendy & Kennicutt|2004^ . In Section [2] we describe the sample. 



observations and data reduction. In Section[3]we discuss the central 
and global stellar populations, and in Section [4] we discuss the ra- 
dial trends and gradients. Full results (i.e. complete radial profiles) 
for the stellar kinematics, gas kinematics, Lick indices, and SSP- 
equivalent population properties are presented in Williams ( 2011) . 



2 OBSERVATIONS AND DATA REDUCTION 

Our sample galaxies are the 28 edge-on galaxies presented in |Bu-| 
|reau & Free man ( 1999) and Chung & Bureau (2004 ). Those works 
describe the sample and observations used here in more detail, but 
in brief half of the galaxies are SOs, half are spirals, and 22/28 host 
a boxy or peanut-shaped bulge. We took long-slit spectra using the 
Double Beam Spectrograph (DBS) on the 2.3 m telescope at Sid- 
ing Springs Observatory. In the instrumental setup used, the DBS 
takes two spectra, each covering a range of ~ 1000 A. The blue arm 
was centred on the continuum absorption features around the Mgb 
triplet. The red arm was centred on the Ha emission line when de- 
tected or moved to the Calcium triplet. The red spectra therefore 



suffer from emission or strong atmospheric absorption and cover a 
wavelength range where stellar population models are poorly con- 
strained. The blue arm spectra are therefore the sole focus of this 
work. 

We reduced the data from the blue arm using standard IRAF 
long-slit techniques, yielding a flat-fielded, wavelength-calibrated, 
sky-subtracted long-slit spectrum along the major axis of each 
galaxy. Cumulative exposure times ranged from 200 to 400 min- 
utes per object. The reduced two-dimensional spectra were linearly 
binned in wavelength and cover the spectral range 4755-57 10 A. 
The spectral resolution of the instrument is 1.2 A full-width half- 
maximum (equivalent to (7j nst « SOkms" 1 at 5150 A). The spatial 
axis has a pixel scale of 0.91 arcsec pixel -1 and covers a maximum 
of 5'. 



Chung & Bureau (2004) have already presented absorption 
line stellar kinematics for these data, but we rederived them for the 
present work for two reasons. Firstly, we need to spatially bin to a 
higher signal-to-noise ratio (S/N) to better measure absorption line 
strengths. Secondly, more powerful and flexible analysis libraries 
and codes are now available, which allow us to treat emission lines 
more precisely (Cappellari & Emsellem 2004, S arzi et aL]|2006) . 
Outside the central few pixels, where the S/N was already large, 
we spatially binned each spectra to a S/N > 30 per A. 

We extracted stellar and gas kinematics using identical tech- 
niques to those described in Williams etaT]{20TTl, i.e. using penal- 
ized pixel fitting (P PXF; |Cappellari & Emsellem 2004 1 and GAN- 
DALF (Sarzi et all|2006) . As in |Williams et al. 1(2011) , we used 
a subset of 88 stars from the MILES library of 985 observed stel- 
lar spectra ( Sanchez-Blazquez et al. 2006b) as templates. We then 
used GANDALF to clean any H/3, [O III] and [N I] emission, which 
are the three main nebular emission features present in our wave- 
length range. The H/3 emission line occurs in the central passband 
of the H/3 absorption Lick index, while [O III] and [N I] occur in the 
continuum passbands of the Fe5015 and Mgb indices, respectively. 
We imposed the kinematics of the [O III] emission line on the H/3 
and [N I] lines. This simplification is justified for our limited goals 
(cleaning emission) because (1) in the high S/N galaxies in our 
sample, where the kinematics of the H/3 emission line can be reli- 
ably constrained, there is no evidence that they differ systematically 
from those of the [O III] line, and (2) there is no evidence that the 
Hj3 and [O III] emission line kinematics differ by an amount greater 
than their uncertainties in similar analyses of either the SAURON 
ETG sample (Sarzi et al.|2006) or the SAURON early-type (Sa/Sb) 
spiral bulge sample {Falcon-Barroso et al.|2006) . 

Using the stellar kinematics and emission-cleaned spectra, we 
use the method of Kuntschner et al. ( 2006 ) to measure the strengths 
of the absorption lines present in our data (H/3, Fe5015, Mgb, 
Fe5270, Fe5335 and Fe5406) in the Lick/IDS system (Burstein| 
[eFaT] [19841 |Worthey| p994[ |Trager et al^[T998l >. We compare 
these Lick index measurements to an interpolated grid of [Thomas^] 
|Maraston, & Bender| (2003) SSP models, yielding SSP-equivalent 
luminosity-weighted ages, metallicities [Z/H] and ce-element en- 
hancement [a/Fe] (using the method of |Proctor et al.|[2004) . We 
exclude the Fe5270 or Fe5335 feature (depending on the galaxy 
redshift) because of bad columns on the CCD of the blue arm of 
the DBS. We exclude Fe5406 in the most distant galaxy (ESO 597- 
G036) because that line is redshifted into the strong [O I] sky line 
at 5577 A. 
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Figure 1. Central H/3, Fe5015 and Mgb Lick indices and luminosity- 
weighted SSP-equivalent population parameters as a function of cen- 
tral stellar velocity dispersion. Indices are measured in magnitudes, i.e. 
— 2.51og| (l — //A/), where / is index value and Al is the width of its 
bandpass, both of which are measured in angstrom. Large colored symbols 
are from our sample; blue squares are our boxy bulges, and red circles our 
round bulges. Error bars are omitted for clarity; the typical uncertainties are 
±0.1 mag in H/3, Fe5015 and Mgb and ±0.1 dex in [Z/H], log(age/Gyr) 
and [a/Fe]. The smaller gray points are comparison data taken from the 
literature. Top two panels: Gray circles are ETGs from Sanchez-Blazquez 
|et~aT1j2006a) and SAURON jKuntschner et al.|2006^ . The shaded green re- 
gions are the regions populated by the bulges of late-type (Sb-Sd) spirals 
in the sample of Ganda et al. 1 2007 1. Bottom three panels: The gray circles 
are ETGs from Thomas et al. 1 2005 1, which, like ours, use the Thomas et al. 
p003) SSP models. 
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Figure 2. Radial [Z/H] gradients as a function of central velocity disper- 
sion. Large symbols are from our sample; blue squares are our boxy bulges, 
and red circles are out round bulges. For clarity, the median uncertainties 
are shown as error bars only on a representative data point. Smaller symbols 
are comparison data; green squares are barred SO-Sb galaxies from |PSB 1 1 1 
gray circles are ETGs from Spolaor et al. (2010 1. The thick lines are the 
mean A[Z/H] for the samples: pale blue for our boxy bulges, pale green for 
fhe |PSBll| barred galaxies, and gray for the |Spolaor et1dT1p)l7)) ETGs. The 
thickness of these lines is is the uncertainty on the mean. 



3 CENTRAL LINE STRENGTHS AND SSP PROPERTIES 

The smallest aperture from which it is meaningful to extract data is 
set by the seeing limit of the observations and the width of the slit, 
i.e. 3" x 1.8". We refer to the average quantities within this aperture 
as 'central'. 

In Fig.[T]we show, as a function of central stellar velocity dis- 
persion (7o, the central values of three representative Lick indices 
and the SSP-equivalent age, [Z/H], and [a/Fe] of our sample galax- 
ies. We also show comparison data for ETGs. We find no evidence 
that the central populations of our barred disk galaxies differ from 
those of ETGs at a given velocity dispersion. 

Almost all of these galaxies are barred, and therefore good 
candidates for observing the effects of secular evolution in their 
stellar populations. At least in terms of their central values, how- 
ever, the stellar populations of these barred SO-Sb disk galaxies 
are no different to those of ETGs (or indeed the six unbarred disk 
galaxies in our own sample). This implies that either secular evolu- 
tion does not affect the stellar populations of the centers of bulges 
of barred SO-Sb galaxies, or its effects are no different to those 
of monolithic collapse and mergers, the putative formation mecha- 
nisms of ETGs. We note that this result is consistent with the find- 
ings of |Thomas~ & Davies (2006 ), who studied the central stellar 
populations of spiral bulges across a broad range of Hubble types, 
but did not consider the role of bars in particular. Since our sample 
is restricted to SO-Sb galaxies, our result does not necessarily con- 
tradict the apparent observation of secular evolution effects on the 
central stellar populations of later-type (Sb-Sd) bulges by Ganda 
[etaLl ( |2007l >. 
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4 RADIAL METALLICITY GRADIENTS 

To measure A[Z/H], we fit a straight line to [Z/H] as a function 
of log/?. We use the full radial extent of the data, which cover, on 
average, the inner ~30 arcsec, i.e. the bulge of these local galaxies. 
We exclude all points inside the seeing limit, which we considered 
separately in the previous section. We plot A[Z/H] as a function 
of central velocity dispersion in Fig. [2] The data for 2 of our 28 
galaxies were only of sufficient S/N to measure the central popula- 
tions, so are excluded from our gradient analysis. For comparison, 
we show the [Z/H] and [a/Fe] gradients from a catalogue of ETGs 
assembled by Spolaor et al. 1 2010). Th e original sources of these 
data are[Proctor|(|2003^,|Brough et al.| p007>, |Reda et al.|p007), 
Sanchez-Blazquez et al. (2007 ), Spola or et al.| 2008 \ , and Spolaor 
et al.| ( |2010) . All are based on long-slit observations, and all sources 
except |Proctor|p003l > use the |Thomas eTaI]p003] l SSP models. We 
note in any case that, while the choice of SSP models may affect the 
absolute stellar population parameters derived, the radial gradients 
are unlikely to be strongly affected ( Kuntsch ner et al.|2010| >. 

From Fig. [2] we conclude the following. (1) There is no evi- 
dence that A[Z/H] in our boxy bulges is correlated with velocity 
dispersion. (2) The boxy bulges of our sample of barred galax- 
ies have shallower metallicity gradients than those of ETGs, both 
on average and at a given velocity dispersion. The mean value of 
A[Z/H] for the boxy and peanut-shaped bulges is —0.06 ± 0.04 and 
there are several cases of positive metallicity gradients. In contrast, 
the mean A[Z/H] of the Spolaor et al. 1 2010) ) catalogue of ETGs 
is —0.23 ± 0.02. These results are qualitatively consistent with the 
simulations of Fried li et al.| (1994) , who find that outflows and 
inflows in barred galaxies make pre-existing radial gradients less 
steep. Since the pre-existing age and [a/Fe] gradients of ETGs and 
bulges are, on average, flat, it is not a surprise that, on average, 
A(age/Gyr) = and A[a/Fe] = for our sample as well. We omit 
these results from Fig. [5] but the results for age and [a/Fe] profiles 
are presented and fully discussed in Williams (20TT). 

One may reasonably worry that line-of-sight effects in our 
sample of edge-on barred galaxies are responsible for some or all 
of the flattening of their radial gradients. While we cannot quan- 
tify this effect, we argue that it must be small for two reasons. 
(1) There is no systematic difference between the radial gradients 
of our 14 SOs (edge-on galaxies largely free of dust) and 14 spi- 
rals (edge-on galaxies with prominent dust lanes). This suggests 
the role of dust is small. (2) As shown in Fig. [2] we see a simi- 
lar result — shallower gradients in barred galaxies — in data taken 
from the bulges of a face-on sample of barred galaxies that cannot 
suffer from line-of-sight flattening (PSB11 1. The statistical signifi- 
cance of the difference between metallicity gradients in the |PSBlT] 
galaxies and unbarred ETGs is admittedly weak (the mean A[Z/H] 
for the |PSB ll] sample is —0.15 ±0.04), but the strong and clear 
correlation between A [Z/H] and a seen in unbarred ETGs with 
log(rj/kms _1 ) < 2.2 |Spolaor et al.|2010} is totally absent from 
both our sample of boxy bulges and the |PSB 1 lj barred galaxies. 



5 DISCUSSION 

In summary, we have shown that the stellar populations at the very 
centres of the bulges found in barred SO-Sb galaxies do not dif- 
fer from those of ETGs of the same velocity dispersion. On larger 
physical scales, however, these bulges do differ from ETGs: they 
lack the correlation between metallicity gradient and velocity dis- 
persion found in ETGs, and their average stellar population metal- 



licity gradients are shallower than unbarred ETGs with the same 
velocity dispersion. 

It is clear that the role of bars in transforming the stellar popu- 
lations of disk galaxies is a significant gap in our understanding of 
galaxy evolution. This work suggests several possible avenues for 
further study. Simulations should attempt to go beyond the quali- 
tative statement that bars should make pre-existing abundance gra- 
dients shallower, either by drawing their initial abundance profiles 
from samples such as |Spolaor et al.| ( |2010| l and Kuntschne r et al.| 
1 2010 1, in which A[Z/H] is correlated with velocity dispersion, or 
by using full cosmological simulations with sufficient resolution to 
capture bar-driven evolution. Simulations should also make predic- 
tions of the vertical gradients of barred galaxies as a function of 
radius (e.g. Friedli 1998 ; William s et al.|20lT) . Integral-field obser- 
vations of the boxy/peanut-shaped bulges of edge-on galaxies are 
the ideal observational tests of such predictions. Observations of 
face-on galaxies are complementary to this work, and integral-field 
data raise the possibility of looking for anticipated azimuthal vari- 
ations in stellar populations, such as local minima in age at the bar 
ends ( |Wozniak|20"07l|Perez et al.|2007| l. 
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